Single-particle spectral function of quarter-filled ladder systems 
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We study the single-particle properties of quarter-filled ladder systems such as a'-NaV20s by 
means of a recently developed generalization of the variational cluster perturbation theory to ex- 
tended Hubbard models. We find a homogeneous antiferromagnetic insulating phase for nearest- 
neighbor repulsions smaller than a critical value, without any metallic phase for small repulsions. 
Different from C-DMFT and LDA considerations, the inclusion of diagonal hopping within a ladder 
has little effect on the bonding bands, while flattening and shifting the antibonding bands. In the 
low-temperature charge-ordered phase, the spectrum depends on whether the ordering is driven by 
the Coulomb repulsion or by the coupling to a static lattice distortion. The small change of the 
experimentally observed gap upon charge ordering implies that the lattice coupling plays an impor- 
tant role in this ordering. Inter-ladder coupling is straightforward to include within our method. 
We show that it has only a minor effect on the spectral function. The numerically calculated spectra 
show good agreement with experimental angle-resolved photo-emission data. 

PACS numbers: 71.10.Fd,71.38.-k 



I. INTRODUCTION 

In recent years low-dimensional strongly-correlated 
systems have been the subject of many experimental 
and theoretical studies due to their fascinating proper- 
ties such as the occurrence of ordered patterns of the 
ion charges. A compound in this class of materials is 
the low-dimensional vanadium bronze a'-NaV2 0s. Al- 
though known for many yearsi it has attracted consider- 
able attention in recent years because of a very interesting 
low-temperature phase. The compound exhibits a spin- 
Peierls-like transition at T c « 35 K accompanied by the 
opening of a spin gapJ> At the same^ or slightly higher 
temperature charge ordering takes place. Different from 
the first X-ray investigations, recent studies^^ at room 
temperature showed a disordered state with equivalent 
valence 4.5 for all vanadium ions. Below the phase tran- 
sition point, NMR-studiesi gave two different valences for 
the ions, a clear evidence for the formation of a charge- 
order pattern below T c . Since one d xy electron is shared 
by two V sites in a V-O-V rung, the ordering occurs 
as a static charge disproportion S between the V ions, 
yielding charges 4.5 ± S with a zig-zag pattern of S's. 
Since the crystal environment of the vanadium ions is 
asymmetric, the d xy electrons are coupled to the lattice 
via a strong Holstein-like electron-phonon interaction: 8 
This results in a static lattice distortion below the charge 
ordering transition temperature, where the ion displace- 
ments from their positions in the high-temperature phase 
is of the order of 0.05 A as observed in X-ray diffraction 
experiments^ 

Although the crystal structure of a'-NaV20s is com- 
posed of nearly decoupled two-dimensional layers that 
consist of coupled two-leg ladders, spin-susceptibility 
measurements* revealed that the system can be rea- 
sonably well described by a one-dimensional Heisenberg 
model. This behavior could be explained by realiz- 
ing that the molecular-orbital state on a rung occu- 



pied by one electron is a key element of the electronic 
structure^ yielding quasi ID magnetic exchange cou- 
plings. In addition angle-resolved photo-emission spec- 
troscopy (ARPES), performed in the disordered high- 
temperature phase showed quasi ID band dispersions of 
the vanadium 3c? bands^ii^ and it was argued that spin- 
charge separation should be present in this system^ 

Previous studiesi^i^iSiii revealed that the electron- 
phonon coupling is very important for the phase transi- 
tion in quarter-filled ladder compounds. For this rea- 
son we study a model Hamiltonian that includes the 
coupling of the d xy electrons to the lattice. The rele- 
vant parameters for the study of lattice effects can be 
obtained from experiments^ (phonon frequencies), and 
from first-principle calculations (lattice force constants 
and electron-phonon coupling) 

Static and dynamic properties of quarter-filled lad- 
der compounds without coupling to the lattice have 
been studied intensively in the past using different 



methods like mean-field approaches 
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exact di- 



agonalization (ED) of small clusters ^^^SL density- 
matrix renormalization-group (DMRG), 27 cluster dy- 
namical mean-field theory (C-DMFT) ^2, and bosoniza- 
tion and renormalization-group techniques^ Recently 
the influence of the lattice coupling on the charge- 
ordering transition was investigated by employing ED 
methods £L 

In this paper we study the single-particle spectral func- 
tion of the compound a'-NaV205, which can be directly 
related to the ARPES experiments, by applying the re- 
cently proposed variational cluster perturbation theory 
(V-CPT)iSI This method is a combination of the clus- 
ter perturbation theory^! and the self-energy functional 
approach (SFA)|22 which provides results for the infi- 
nite lattice and allows to study symmetry- broken phases. 
It was used with success for the investigation of the 
magnetic ground-state properties of the two-dimensional 
Hubbard models For Hamiltonians including offsite 



Coulomb interactions an extension of this theory has 
been developed^ which turned out to give very accu- 
rate results for the one and two-dimensional extended 
Hubbard model. 

The paper is organized as follows. In Sec. ^ we intro- 
duce the model Hamiltonian and give a short description 
of the V-CPT method. Sec.lrHland Sec .[El includes our 
results for single and coupled ladders, respectively, and 
we finally draw our conclusions in Sec.lVl 

II. MODEL AND METHOD 

On a microscopic level, a'-NaV205 can be described by 
an extended Hubbard model (EHM). In order to take into 
account lattice effects we further extend this well-known 
model by a Holstein-like electron-phonon coupling, which 
a recent LDA studyi 9 - showed to be especially important, 
yielding the model 

H = -Hehm + -Hi + -He-i) (1) 

with i?EHM the EHM Hamiltonian, H\ the contribution 
of the lattice, and H e -\ the Holstein coupling. These 
terms are given by 

#EHM = - **i ( C ia C j<7 + Rc -) 

+ U y^n^riii + 2J VijUiUj, (2a) 

ffi=«E| ( 2b ) 

i 

H e -i = -Cj2^i, (2c) 

i 

where (ij) connects nearest-neighbor bonds, and Uj is 
the corresponding hopping matrix element. In Fig.^ 
the lattice structure and the hopping processes used in 
this study are shown. The most commonly used set for 
these matrix elements is t a = 0.38 eV, % = 0.18 eV, and 
t X y = 0.012 eV, and was obtained by fitting the LDA 
bands£ A recent studj^ gave similar parameter values. 
By including the additional hopping term td in a massive 
downfolding procedure, Mazurenko et. found simi- 
lar values for t a and t xy , but the values t& = 0.084 and 
td = 0.083 differ considerably from previous studies. In 
the present study we set t a as the energy unit and fix 
tf, = 0.5, except for Sec. lIII 01 where we study the spec- 
tral function including the hopping term td- The on- 
site Coulomb interaction is set to U — 8 throughout the 
paper, in accordance with band-structure calculations, 19 
and the inter-site Coulomb interaction Vy is treated as a 
free parameter of the system, since the determination of 
a proper value within first-principle calculations is very 
difficult. 

The lattice distortions are given in units of 0.05 A, 
since the ion displacements in the ordered phase are of 
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FIG. 1: Clusters used for the V-CPT calculations. The diag- 
onal hopping td is indicated only once, but is present equiv- 
alently between other sites. Decoupled bonds treated pertur- 
batively are marked by dashed lines, and the boxes show the 
clusters of finite size. Left: Single ladder with 6x2 cluster. 
Right: Super cluster consisting of two 12 site clusters. 

this order of magnitude. The electron-phonon coupling 
C and the lattice rigidity k were determined by first- 
principle calculations^ 9 - yielding C — 0.35 and K — 0.125 
in these units. We restrict our investigations to static 
distortions and neg lect dynamical phonon effects, similar 
to parts of Ref.[I3. Moreover, as discussed below, we use 
a staggered zig-zag configuration for the Zj, as observed 
experimentally. 9 

The method we use in this paper for the calculation 
of the single-particle spectrum is the variational cluster 
perturbation theory for an EHMi^i The main idea is to 
decouple the lattice into clusters of finite size as depicted 
in Fig.^ yielding the Hamiltonian 

H = J2[ H o C \R) + Hi(R)} + £ H®(R,R'), (3) 

R R,R' 

where R denotes the individual clusters. The first sum 
consists of decoupled intra-cluster Hamiltonians with in- 
teraction part Hi(R), and the second sum gives the cou- 
pling between clusters, which must be of single-particle 
type and is of the general form 

H^(R,R')=J2 T T' c ka^, b - W 

a . b 

The indices a and b are orbital indices within a cluster. 

In the case of the EHM this decoupling into clusters 
yielding Eq. @ cannot be done in a straightforward way, 
since the Coulomb interaction on the decoupled bonds is 
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of two-particle type. In order to get a Hamiltonian of the 
form Eq. © it is necessary to do a mean-field decoupling 
of the Coulomb-interaction terms on inter-cluster bonds, 
which is described in detail in Ref.0- The coupling of 
the clusters is then again of single-particle type, Eq. J3J, 
but in addition this mean-field approximation introduces 
onsite potentials Aj on the cluster boundary, which cor- 
respond to the average electronic densities on site i of 
the cluster boundary. Let us stress at this point that the 
\ are not variational parameters in the sense of V-CPT 
but external parameters to the Hamiltonian Eq. in its 

(c) 

mean-field decoupled form, entering term Hq . 

After the mean-field decoupling of inter-cluster bonds, 
we can apply the V-CPT to the Hamiltonian Eq. ©. Ob- 
viously the Hamiltonian is invariant under the transfor- 
mation 

4 c \R)^Hf\R)+0(R) 
Hp(R,R') - Hp(R,R')-6 R , R >0(R), 
with an arbitrary single-particle operator 

0(fl)=5>a,> c k- C *.*- (6) 

a , 6 

This transformation allows the study of symmetry- 
broken phases by the inclusion of fictitious variational 
symmetry-breaking fields, which do not change the over- 
all Hamiltonian, Eq. but instead just rearrange it. 
They are therefore different in character from the exter- 
nal lattice distortions, Eq. J2bj), and also from the exter- 
nal mean-field parameters A;, which directly enter the 
Hamiltonian. For our model, symmetry-breaking fields 
corresponding to charge order and to magnetic order 
can be important. We will later see that we can omit 
a charge-order symmetry-breaking field. A symmetry- 
breaking field for magnetic order will be considered in 

sec. inrei 

After introducing the external mean-field parameters 
Xi and the variational parameters A = A a .h, they have 
to be determined in a proper way, which is done within 
the framework of the SEA. It provides a unique way to 
calculate the grand potential of a system by using dynam- 
ical information of an exactly solvable reference system, 
which is in our case the decoupled cluster. This grand 
potential is parameterized as a function of the external 
mean-field parameters Xi and the variational parameters 
A = A a The f unctio nal form of the grand potential is 
taken from Refs. I32l34i 

n(A,x l ) = n'(A,x l ) 

+ T V tr In — — 

G<°J(A ij ««; n )- 1 -E(A 1 A< > ta;„) 

-LT^trln(-G'(A,A„i W „)), (7) 

where fV(A, Xi) is the grand potential of the de- 
coupled clusters, G^ (Xi, iui n ) is the non-interacting 



Green's function of the original infinite-lattice prob- 
lem after mean-field decoupling of inter-cluster Coulomb 
interactions, H(A,Xi,i(J n ) is the cluster self energy, 
G'(A, Aj, ioJn) the cluster Green's function and L denotes 
the number of clusters. All cluster properties can easily 
be calculated by the Lanczos algorithm. The sum over 
Matsubara frequencies in Eq. JjJ is evaluated by a contin- 
uation to the real frequency axis, iuj n — > w + iO + , yielding 
an integral from minus infinity to the chemical potential 
fi, determined below. Note that the Hamiltonian, Eq. QJ, 
does not involve /x in any way. 

The general procedure to determine the Aj and the 
values of the variational parameters A a ^ is the follow- 
ing. First one has to distinguish between the external 
parameters Xi and the variational parameters A a) (,. For 
the latter ones the general variational principle of the 
SFA says that f2, Eq. J7J, must have a stationary point 
with respect to A a j, but the SFA does not provide any 
information on the second derivative. That means that 
this stationary point can be a maximum, minimum, or a 
saddle point. The situation is different for the external 
parameters. Since the Aj are mean-field parameters, one 
has to look for a minimum of f2 with respect to the Aj. 
Finding the minimum in the grand potential is equiva- 
lent to a self-consistent solution for the Xi, as shown in 
the appendix of Ref.l34l 

In practice we used the following procedure. For a 
given value of the Xi one has to find the stationary point 
of fi(A, Aj) with respect to A, yielding a function = 
fi(Aj). The proper value of the external parameters Aj is 
then given by the minimum of this function. 

In Sec. IIII Dl we study the effect of the lattice dis- 
tortions Zi, and the optimal distortions are determined 
by the minimum ground-state energy. From a technical 
point of view these distortions are treated on the same 
level as the mean- field parameters Aj, since they are ex- 
ternal parameters to the Hamiltonian, Eq. Q, as well. In 
this case on has a function fi(A, Aj, zi). Again, for each 
pair Aj , Zj one looks for the stationary point with respect 
to A, and the proper choice of Aj and Zi is then given by 
the minimum of the function f2(Aj, zi). 

The single-particle Green's function is then calculated 

by 

G,( W )=[G'H- 1 -T,]- 1 (8) 

with the Fourier-transformed matrix elements T q a 
After applying a residual Fourier transformational one 
finally obtains the fully momentum-dependent Green's 
function G(k,ui) for the infinite size system. Note that 
the external parameters (Xi,zi) are only present in the 
calculation of G (cj), whereas the variational parameters 
A also enter T q , see Eqs. (|4I5I6I) . 

Since calculations are not done at half filling, the chem- 
ical potential /i is not known a priori. However, the 
knowledge of /i is important for the evaluation of the 
grand potential, as discussed above. One can calculate fi 
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from the condition 



It. 



dtoA(k, u>), 



where the spectral function A{k,uS) is given by 



A(k, u) = Im G(k, u> + irj) 

7T 



(9) 



(10) 



and r\ is a small Lorentzian broadening. This amounts 
to a self-consistent procedure, since for the calculation 
of fi the Green's function G(k,u>) is needed, and for the 
determination of G(k,uj) one has to know ji. This cycle 
can be avoided as follows. One can infer the chemical 
potential directly from the energies of the excited states 
obtained by the ED. An approximate value for the chem- 
ical potential is given by 



PIPES I pPES 



MED 



(11) 



with E^^ s the minimal energy of inverse-photo-emission 
(IPES) states and the maximum energy of photo- 

emission (PES) states. This value only weakly depends 
on the mean- field and variational parameters. As dis- 
cussed below, in all our calculations we found a well es- 
tablished gap between the PES and IPES states yield- 
ing a constant density n in a reasonably large neighbor- 
hood of the physical chemical potential /i, in agreement 
with quantum Monte Carlo calculations^ Therefore /ied 
gives a reasonable approximation for our calculations. 

For the details of this method and the calcu lation of the 
grand potential we refer the reader to Refs. I3fl3 1I33I34I 
and references therein. 



III. RESULTS FOR SINGLE LADDERS 

A. Critical coupling 

We start our investigations with decoupled ladders, 
i.e., t xy = and V xy = (Fig.[TJ. Before we turn to 
the spectral function, we study the charge-ordering tran- 
sition as a function of V = V a = Vb- Let us first look 
on the effect of the mean-field decoupling of the non- 
local Coulomb interactions across cluster boundaries as 
described in Sec.[n] For this purpose we consider the 
EHM without coupling to the lattice in the limit of ex- 
actly one electron per rung and tb = 0, i.e. without hop- 
ping between rungs. In mean-field approximation this 
case results in a second-order phase transition between 
a disordered state and a zig-zag ordered state at a crit- 
ical interaction of V C MF = 1.0. On the other hand this 
case is exactly solvable by a mapping to an Ising model 
in a transverse field^ yielding a critical interaction of 
■j/exact _ 2.0i2i Thus we expect strong mean-field ef- 
fects, since in this special limit we found V^ &ct = 2V C MF . 
Since it can be assumed that a finite value of the hopping 



between adjacent rungs tb weakens the charge ordering, 
the actual critical value V c is presumably located slightly 
above 2.0 when tb is included. 

For single ladders at quarter filling one has only two 
different values for the mean-field parameters A, , namely 
Xa on sublattice A and Xb on sublattice B, which corre- 
spond to the left and right side of the rung on the cluster 
boundary, respectively, see Fig.^ In order to reduce the 
number of mean-field parameters, we set Xa = {n}+6 and 
Xb = i n ) — S, with the average density fixed to (n) = 0.5. 
This yields only one mean-field parameter S instead of 
two parameters, A ,4 and Xb- 

For second-order phase transitions it can in addition 
be important to rearrange the Hamiltonian by means 
of a fictitious staggered chemical potential as a varia- 
tional parameter t24 This field is included by adding and 
substracting (Eq. (0) the single-particle operator O(R), 
Eq. |JBJ, with 



iQr a 



(12) 



where e denotes the variational parameter, r a is the lat- 



tice vector of site 



is the Kronecker-<5 for sites a, b, 



and Q = (tt, it). Initial calculations showed that in the 
present case the inclusion of such a field does not have 
any significant effect and the relative change in f2 is at 
most of the order of 10 -4 . Although the staggered field 
Eq. ljT2*|) is the Weiss field associated with the CDW, its 
effect is very small, since the symmetry is already bro- 
ken by the mean-field decoupling at the boundaries. We 
also considered a staggered chemical potential similar to 
Eq. I|12l) , but which is nonzero only on the cluster bound- 
aries. In the spirit of V-CPT this corresponds to include 
the mean-field parameter 5 also in the set of variational 
parameters A. Not surprisingly, the effect of this field 
was even smaller than for the field Eq. (|12f> . For this 
reason all further calculations have been done without a 
staggered chemical potential. 

We now discuss our results, setting tb — 0.5. We 
first discuss the case without coupling to the lattice. 
Finite lattice coupling will be considered in Sec. lIII Dl 
In order to determine the order of the transition within 
the framework of V-CPT it is sufficient to calculate the 
grand potential f2, Eq. 0, as a function of the mean- 
field parameters^ 4 - Fig.0 shows the dependence of f2(<5) 
on the mean-field parameter 6 calculated with a 6 x 2 
cluster as reference system. One can see that the system 
undergoes a continuous phase transition^ which is lo- 
cated between V = 1.5 and V — 1.7, since the minimum 
of Q(S) shifts from S = to a finite value. This value for 
the critical interaction is considerably smaller than the 
above mentioned value of the analytical solution, but the 
agreement is much better than the result of the purely 
mean-field calculation, V^. MF = 1.0. 

In order to study the finite-size dependence of the 
critical Coulomb interaction we performed calculations 
on clusters of different length, and the results are de- 
picted in Fig.|3| The steps in V in our calculations were 
AV = 0.01, which results in error bars of AVv = 0.005. 
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FIG. 3: Finite-size dependence of the critical Coulomb inter- 
action V c without lattice coupling. Error bars are due to the 
finite step AV = 0.01 in the calculations. Dotted line: Lin- 
ear extrapolation of the 8 and 6 rung cluster. Dashed line: 
Quadratic extrapolation of the 8, 6, and 4 rung cluster. 



As expected, V c is strongly finite-size dependent. From 
Fig. |3 we can expect that for larger cluster sizes the crit- 
ical interaction V c increases further and reaches the ex- 
pected value of slightly above 2.0, but for a more sophis- 
ticated finite-size scaling our cluster sizes are too small. 
Nevertheless it is possible to study the spectral function 
both in the disordered and the ordered phase. Since the 
calculations for the 8x2 ladder are very time consum- 
ing, all single-ladder spectra presented in this paper have 
been determined with a 6 x 2 ladder as reference system. 



FIG. 4: Single-particle spectral function A(k,u>) calculated 
with a 6 x 2 cluster in the disordered phase at V a — 14 = 1.3. 
Top panel: Momentum k a = perpendicular to the ladder. 
Bottom panel: k a = 7T. The dashed line marks the chemical 
potential calculated by Eq. (f 1 1 pi , the dotted line marks the 
result obtained from Eq. 



B. Disordered phase 

We start our investigations of the spectral function 
with the disordered high-temperature phase. Since a!- 
NaV205 may be near a quantum critical point between 
ordered and disordered phase, we choose the nearest 
neighbor interaction to be slightly below the critical 
value. We set V = V a = Vb = 1.3, tb = 0.5, and we do not 
include diagonal hopping, i.e. td = 0. The result of this 
calculation is shown in Fig.^ An additional Lorentzian 
broadening of 77 = 0.1 has been used for all spectra shown 
in this paper. The dashed vertical line marks ^ed cal- 
culated from Eq. ((Til , and the dotted line denotes ^cpt 
determined from the condition Eq. @ ■ For the latter 
quantity the sum over momentum vectors had to consist 
of about 80 vectors in order to get a well converged re- 
sult. It is easy to see that /xed = 1.71 lies exactly in 
the middle of the gap, whereas ^cpt — 1-23 is located at 
its lower boundary. But since there are no in-gap states 
both values of \x give approximately the same average 
density n, and the ground-state energy Eq = f2 + fJ-N 
hardly depends on whether we use /j,et> ot Mcpt- These 
facts confirm that our approximation to use /x = ^ed as 
chemical potential gives correct results, and in addition 
the numerical effort for this procedure is much less than 
for the above described self-consistent determination of 
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FIG. 5: Gap A in the spectral function as a function of V a - 
Squares: Vb = 0. Diamonds: Vb = V a . 

As one can easily see in Fig.0] the spectral function 
exhibits a well defined gap around the chemical poten- 
tial, a clear indication of insulating behavior. In order 
to check if the insulator is only stable above some criti- 
cal inter-site Coulomb interaction, we calculated the gap 
A at (k a ,kb) — (0, 7r/2) as a function of the intra-rung 
interaction V a . We studied two cases with Vb = and 
Vb = V a , respectively, and the results are shown in Fig.0 
Note that for Vb — no mean-field decoupling is needed, 
since there are no interaction bonds between different 
clusters. At V a = 0, where both cases are equivalent, we 
found a finite value of the gap, A w 0.55. We checked the 
finite size dependence by calculating the gap on a 4 x 2 
cluster giving A as 0.59. By applying a linear l/iV run g S 
extrapolation to N rungs = oo one gets A rj 0.47, indicat- 
ing that the curves in Fig. [S] somewhat overestimates the 
value of the gap for the infinite ladder. Nevertheless we 
conclude from our calculations that the system is insulat- 
ing already for small values of V a ■ This is consistent with 
DMRG calculations^— where for t a > % a homogeneous 
insulating phase has been found for V = 0. The behav- 
ior of the spectral function is also in agreement with ED 
calculations on small clusters for V = 0, where for large 
enough t a an insulating state has been founds Similar 
results have been 

obtained by Kohno^i for the U = oo Hubbard ladder. 

In the case Vb = 0, which means that there is no 
Coulomb interaction between adjacent rungs, we found 
that A increases linearly with V a ■ For Vb = V a the gap is 
slightly larger and the deviation increases with increas- 
ing V a - Here, with a 6 x 2 cluster as reference system, 
the system starts to order at V c ks 1.625, which results 
in the kink in A around this critical value. Note that for 
Vb = such a phase transition is not possible. 

Let us now discuss the spectral features for k a = as 
shown in Fig.^J The spectral function looks very sim- 
ilar to that of the half-filled one-dimensional Hubbard 
model with a totally filled lower and an empty upper 
band. Different from the ID Hubbard model the gap 




h 

FIG. 6: Magnetic properties of a single ladder at V a = Vb = 
1.3. Upper panel: Spin correlation function S r = (SiSi +r ) 
calculated on an isolated 6x2 cluster. Lower panel: Grand 
potential fi as a function of the strength of the fictitious field 
Ea.lTSl. 



between these two bands is not only determined by the 
onsite interaction U, but mainly by the intrarung inter- 
action V a , as discussed above. At k = (0, 0) one can see 
signatures of spin-charge separation, where the band is 
split into a low energy spinon band (at approximately 
u> — fi ~ — 1.5) and a holon band at slightly higher en- 
ergy (u> — fx « —2.0). This splitting has not been seen 
directly in experiments^*^ since it is small and tem- 
perature effects did not allow a high enough experimen- 
tal resolution. However, by studying the temperature 
dependence of ARPES spectra, it was argued that sub- 
tle spectral-weight redistributions can be related to spin- 
charge separation^ Some spectral weight can also be 
found at very high energies of about (u> — fj,) « 8.5, which 
is close to the onsite energy U — 8 and can thus be re- 
lated to doubly-occupied sites. 

Infrared (IR) experiments probe transitions near the T 
point, that is between even (0,0) and odd (0,7r) states 
in the language of single ladders. From Fig.0] one can 
extract an excitation energy of roughly 3i a , which is in 
good agreement with the experimentally found 1 eV ab- 
sorption peakii^ 

The lower and upper band disperse with period 7r in- 
dicating a doubling of the unit cell in real space, sim- 
ilar to the ID Hubbard model. In order to determine 
the origin of this doubling we have calculated the real- 
space spin correlation function S r — (SfSf +r ) within the 
cluster by exact diagonalization, where S* and S* are 
the ^-components of a spin on the cluster boundary and 
on a rung with distance r to the boundary. In the up- 
per panel of Fig.Elthis correlation function is shown for 
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two different paths, where the solid line is S r along one 
leg of the ladder, and the dashed line is S r on a zig- 
zag path through the ladder. Both correlation functions 
show clear antiferromagnetic correlations along the lad- 
der similar to results obtained by the finite-temperature 
Lanczos method. 25 By applying a fictitious symmetry- 
breaking magnetic field via the operator O(R), Eq. JfJJl, 
we can estimate whether this ordering is of long-range 
type or not. Similar to Ref.l33l we choose for this field 

A a , 6 = hS a , bZ(T e^ ra , (13) 

where z a is ±1 for spin T, J. in the orbitals a,b, 8 a . b is 
the Kronecker-<5, and h is the field strength. The wave 
vector Q is set to (0,tt) yielding a staggered field along 
the ladder. The dependence of f2 on this fictitious field 
is depicted in the lower panel of Fig-EI Similar to the 
one-dimensional Hubbard model at half filling^ there is 
only one stationary point at h = 0, which means that 
the system does not show long-range antiferromagnetic 
order, but is rather in a paramagnetic state with short- 
range antiferromagnetic correlations. 

The above considerations show that the system ex- 
hibits short-range antiferromagnetic spin correlations 
along the ladder, which can produce the doubled unit 
cell. Nevertheless it is also possible that the doubling 
of the unit cell is due to short-range charge correlations 
and not due to spin correlations. In order to clarify this 
point we calculated the spectral function at Vb = and 
finite V a , where no charge ordering is possible. Also in 
this case the periodicity of the bands with largest spec- 
tral weight is tt at k a — and 2tt at k a = tt. This shows 
that below the phase transition at V a = 1.3, the doubling 
of the unit cell is mainly due to short-range spin corre- 
lations, and charge correlations play only a minor role in 
this context. 

When turning to k a = tt, the spectral function looks 
totally different. As on can easily see in Fig.0] there is 
hardly any spectral weight below the chemical potential, 
which means that there are no occupied states in the 
channel k a = tt. This can be understood, because k a = 
tt corresponds to an antibonding state within a rung, 
which has energy 2t a relative to the bonding orbital and 
is therefore not populated in the ground state. 

An obvious difference between the spectra for k a = 
and k a — tt is that in the latter case the excitations 
with largest spectral weight located between uj « 3 and 
lu w 4.5 disperse with periodicity 2tt instead of tt. Quali- 
tatively this can be understood as follows. When insert- 
ing a particle with k a = 0, this electron will occupy a 
state in the bonding orbital. Since one of the two states 
in this orbital is already occupied, the additional par- 
ticle must have opposite spin, and thus this particle is 
connected to the antiferromagnetic background. A par- 
ticle with k a — tt occupies a state in the antibonding 
orbital, and since this orbital is not occupied, both spin 
directions possess equal possibility. Therefore an electron 
with k a — tt is not influenced by the antiferromagnetism 
in the ground state. 
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FIG. 7: Spectral function A(k,ui) when the diagonal hopping 
is included, tt = 0.25, td = 0.25. The Coulomb interaction 
was V a = Vb = 1.3. The dotted line marks the chemical 
potential. 



C. Disordered phase including diagonal hopping 

So far we have studied single ladders only with hop- 
ping parameters t a and th and neglected additional diag- 
onal hopping processes td, as indicated in Fig.^ These 
hopping processes have been important in first-principle 
calculations in order to fit the LDA bands correctly^ 
Moreover td has been important in C-DMFT calculations 
in order to describe the insulating state in the disordered 
phase properly^ In this section we study the effect of td 
within the V-CPT framework. 

In Fig.[7|the spectral function is shown for V a — V = 
1.3, tb = 0.25 and td = 0.25, where the hopping parame- 
ters are chosen similar to Refill Whereas the spectrum 
at k a = is almost indistinguishable from Fig. 01 we see 
a big difference at k a — tt. There is still hardly any spec- 
tral weight below the Fermi energy but the band with 
largest spectral weight above the Fermi level is now lo- 
cated at approximately w — [i w 2.0 and can be regarded 
as dispersionless. 

^From a qualitative point of view this can be explained 
by the dispersion of non-interacting fermions on a two- 
leg ladder in the presence of diagonal hopping, which is 
given by 

e[k) = - t a cos k a - 2t b cosk b 

(14) 

— 2^cos k a cos k b , 

where the values for k a are restricted to and tt, and in 
these two cases the dispersion can be written explicitly 
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as 

e(k a = 0, k b ) = -t a -2(tb + t d ) cos k b (15a) 
e(k a = 7r, kb) = +t a -2(t b - t d ) cos k b . (15b) 

This means that for k a = the bandwidth is determined 
by the sum of t b and td, whereas for k a = tt it is set by the 
difference of these two hopping processes. Since we used 
tb = td = 0.25, this fits perfectly to the spectrum shown 
in Fig.0 The sum is equal to the value of t b used for 
the calculations without diagonal hopping in Sec. 1111 Bl 
and the difference is equal to zero, which explains the 
dispersionless band at k a = tt. 

The picture that evolves from our calculations is some- 
what different to that obtained in first-principle and C- 
DMFT calculations. To begin with, the bands obtained 
from the LDA all disperse with periodicity 2ir and not it, 
as observed experimentally, along the b direction. More- 
over we could not find any signature of a flattening of 
the upper d xy bands in the direction k — (0, 0) — > (0, tt) 
when a diagonal hopping is included, which was reported 
in Ref. |33- The main difference of our calculations to C- 
DMFT results is that C-DMFT finds a metal-insulator 
transition at some finite value of V, and this transi- 
tion point is shifted downward significantly when td is 
included^ In contrast we find an insulating state at rea- 
sonable values of U already for V = 0, even without the 
inclusion of td- The discrepancy to C-DMFT calculations 
are very likely due to the fact, that the cluster used in 
the C-DMFT calculations consisted only of a single rung, 
and fluctuations along the ladders, which seem to be im- 
portant in this system, have been neglected altogether. 



D. Ordered phase 

We investigate two different driving forces for the oc- 
currence of a charge-order pattern, (i) the coupling of the 
electrons to lattice degrees of freedom, and (ii) nearest- 
neighbor Coulomb interaction, similar to Ref.ll7l 

Let us start our investigations with possibility (i), the 
coupling to the lattice. In order to keep the calcula- 
tions simple, we consider static lattice distortions, as dis- 
cussed in Sec.[n] The inclusion of dynamical phonon ef- 
fects would pose a severe problem to the diagonalization 
procedures, because for phonons the Hilbert space is a 
priori of infinite size, and some truncation scheme has 
to be appliediii Well converged results for the spectral 
function in the presence of dynamical phonons have so 
far only be achieved for the polaron and the bipolaron 
problem^ 

We assumed a zig-zag charge order pattern, justified 
by experimental evidenced 

z t = ze lQr '. (16) 

In order to keep the calculations simple and the num- 
ber of independent variables small, we did not consider an 



additional variational parameter like a staggered chem- 
ical potential. The proper value of the distortion z is 
determined as discussed in SeclLTl 

Motivated by previous workjil we set V a = V b = 1.3, 
since for this choice we expect the distortions to be close 
to the experimental value of z cxp ~ 0.95 (in units of 
0.05 A). Indeed we found z — 0.911, which is close to 
z OX p, and the mean-field parameter was <5 = 0.338. The 
spectrum calculated with these values is shown in the top 
plots of Fig.GD 

The spectral function shows similar features as in the 
undistortcd phase. For k a = the bands disperse with 
periodicity tt, whereas for k a — tt no evidence for a dou- 
bling of the unit cell can be found, and the periodicity is 
2ir. Nevertheless, the gap at k = (0, n/2) is considerably 
larger than for V = 1.3 without distortions, see FigJIj 

An interesting quantity when considering charge- 
ordering phenomena is the charge order parameter, which 
we calculate as 

meow = j^rr £ «n 3 -> - <n» , (17) 
j 

where the expectation value (rij) is calculated from the 
V-CPT Green's function, and (n) = 0.5. The factor (n) 
in the denominator assures that the order parameter is 
normalized to the interval [0,1]. For V a = Vb = 1.3 
and static distortions we obtained tocdw = 0.65, which 
means that the disproportion of charges is rather large. 

Let us now consider possibility (ii) , where the coupling 
to the lattice is switched off, C — 0, n = 0, and the charge 
ordering is driven by the nearest-neighbor Coulomb inter- 
action. In order to make a connection to the results ob- 
tained with lattice distortions, we calculate the spectral 
function at a similar value of the charge order parameter. 
We found that for V a — V b = 2.05 the order parameter is 
tocdw = 0.66, close to the value found above. 

The spectral function is shown in the lower plots of 
Fig-IHl The spectral features again look very similar to 
Fig.|H By comparing the upper and lower plots of Fig.|Hl 
one can see that in both cases the gap at k = (0, tt/2) is 
larger than in the disordered phase, Fig.Q] To be specific 
we found a gap size of approximately 2.6t a in the presence 
of lattice distortions and 3At a without lattice distortions, 
whereas in the disordered phase the gap was l.8t a - It is 
interesting that the momentum-resolved single-particle 
spectral features do depend on the driving force of the 
transition, which was much less pronounced for, e.g., spin 
and charge susceptibilities obtained by integration over 
the electron states 

The excitation energy near the T point, relevant for 
IR experiments, can be read off from Fig.[S]to be roughly 
At a with lattice distortions and 5t a without distortions. 
Although these excitation energies are not constant com- 
pared to the disordered phase, calculations including the 
lattice degrees of freedom give a better agreement to ex- 
perimental IR absorption datajS which show neither a 
shift of the 1 eV peak nor the appearance of new peaks 
related to electronic transitions. 
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FIG. 8: Spectral function A(k, u) in the ordered phase. Top: Transition driven by coupling to the lattice. Bottom: Transition 
driven by Coulomb interaction. For the interaction V = V a — Vb was used, with the values as given in the plots. The dotted 
line marks the chemical potential. 



IV. RESULTS FOR COUPLED LADDERS 

So far we studied only single ladders and neglected 
the inter-ladder couplings, since they are frustrated and 
one might assume that they are only of minor impor- 
tance. Nevertheless our approach allows us to include 
these inter-ladder couplings by choosing an appropriate 
cluster geometry, as indicated on the right hand side of 
Fig.Q] Note that it is necessary to use a 2 x 12 super- 
cluster which allows for a commensurate charge order 
pattern across the cluster boundaries. For details of the 
treatment of super-clusters we refer the reader to Ref.l34l 

The parameter values for the inter-cluster coupling are 
chosen in the following way. First-principle calculations 
have shown that the effective hopping between different 
ladders is very small, so we set t xy — 0.1t a , and longer 
ranging hopping processes are neglected since the linear 
dimensions of the cluster are rather small. The values 
for the other parameters are the same as used for the 
calculations in Sec. 1111 fcfl 

In Fig.|51the spectral function for V a = Vb = V xy = 1.3 
is shown. For k parallel to the b axis one can easily see 
that the spectrum looks very similar to the spectrum of a 
single ladder (upper panel of Fig.^J . The main difference 
between single and coupled ladders is that the chemical 
potential is much larger in the latter case (/i w 3.0), 
which is due to the frustrated inter-ladder bonds. 

When turning to k parallel to the a axis the spec- 
tral function looks very different. The most striking fea- 



ture is that there is hardly any dispersion of the bands, 
and the filled low-energy band can actually be considered 
as dispersionless. The spectral weight of this excitation, 
however, decreases significantly away from k = and is 
transfcrcd to unoccupied states above the Fermi level at 
approximately u> — ji i=s 1.5. 



Let us now compare our numerical results to experi- 
mental data. Kobayashi et. a! 11 performed angle-resolved 
photo-emission spectroscopy (ARPES) at room temper- 
ature, where the system is in the disordered phase. For 
momentum transfer parallel to the a direction they found 
no dispersion of the V 3c? bands, which fits perfectly well 
to our results. For k along the b axis a band dispersion of 
a ID antiferromagnetic quantum system was found with 
experimental band with of approximately 0.06-0.12 eV. 
This value is rather small compared to the band width 
in our calculation of approximately 0.35 eV, see Fig. [5] 
We checked that the band width scales with the hop- 
ping along the ladder tb (not shown), and therefore this 
discrepancy between calculations and experiment can be 
shortened by choosing a smaller value for tb, which does 
not significantly affect the charge ordering of the system. 
Nevertheless the strong difference between spectra along 
a and b direction are well described by our calculations. 
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FIG. 9: Spectral function A(k,ui) in the disordered phase at 
V a = Vb = V xy — 1.3 calculated on the 2 x 12 super-cluster. 
Top: Momentum k along the ladder direction. Bottom: k 
perpendicular to the ladder direction. 



V. CONCLUSIONS 

In this paper we have applied the recently proposed 
generalization of the V-CPT for extended Hubbard mod- 
els to the case of quarter-filled ladder compounds. We 
were thus able to perform the first theoretical study of 
the spectral function of a'-NaV205 within the extended 
Hubbard model. 

For single ladders in the disordered phase we found 
that in the channel k a = the system behaves like a 
one-dimensional antiferromagnetic insulator, and the gap 
is mainly determined by the nearest-neighbor Coulomb 
interaction on a rung. Our calculations suggest that the 
system is in an insulating phase for all values of V . 



This picture still holds when a diagonal hopping td 
is included in the Hamiltonian, which was suggested to 
be important by LDA and C-DMFT studies. We could 
show that for k a = hardly any changes can be seen 
in the spectral function, whereas for k a — tt the bands 
become flat. These findings do not agree with LDA 
considerations^ where a flat upper band was observed 
for k a = 0, thus requiring a finite value of td in a tight- 
binding fit. 

For the transition into the charge-ordered low- 
temperature phase we considered two driving mecha- 
nisms, the coupling to a static lattice distortion and the 
nearest neighbor Coulomb interaction. With lattice cou- 
pling we found, similar to Ref.IrR that for V — 1.3 the 
distortion is close to the experimentally found size, ac- 
companied by a large disproportion of charges. In or- 
der to reach the charge-ordered phase solely by Coulomb 
interactions, we had to use a large value of V (with 
V = 2.05 for the same value of the order parameter), 
which resulted in a large gap in the spectral function, 
considerably larger than in the disordered phase and 
in the ordered phase with lattice distortions. Since IR 
experiments^ do not show such a discrepancy, we sug- 
gest that for the description of the ordered phase, lattice 
distortions cannot be neglected. 

Within our approach it was straightforward to study 
the effects of inter-ladder coupling on the spectral func- 
tion. We found that the spectra along the ladder di- 
rection are not significantly affected by these couplings. 
Perpendicular to the ladders the calculated bands are al- 
most dispersionless, in good agreement with experimen- 
tal data. 
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